On continuum modeling of sputter erosion under normal incidence: 
interplay between nonlocality and nonlinearity 
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Under specific experimental circumstances, sputter erosion on semiconductor materials exhibits 
highly ordered hexagonal dot-like nanostructures. In a recent attempt to theoretically understand 
this pattern forming process, Facsko et al. [Phys. Rev. B 69, 153412 (2004)] suggested a nonlocal, 
damped Kuramoto-Sivashinsky equation as a potential candidate for an adequate continuum model 
of this self-organizing process. In this study we theoretically investigate this proposal by (i) formally 
deriving such a nonlocal equation as minimal model from balance considerations, (ii) showing that 
it can be exactly mapped to a local, damped Kuramoto-Sivashinsky equation, and (iii) inspecting 
the consequences of the resulting non-stationary erosion dynamics. 

PACS numbers: 68.55.-a,79.20.-m, 02.60.Lj 



I. INTRODUCTION 



This equation is given by 
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Sputter erosioni, the bombardment of solid target sur- 
faces with ionized particles to remove or to detach target 
material, has a long tradition in physics as an experimen- 
tal technique to clean, smooth, or appropriately prepare 
solid surfaces. The preparation of on the nanoscale per- 
fectly flat surfaces, however, does not seem to be possible. 
This is because various surface roughening and smooth- 
ing processes compete during the sputtering process and 
can lead to self-organized pattern formation of the eroded 
surface morphology. Depending on the target material, 
temperature, ion beam energy, angle of incidence and 
various other parameters, one generically observes the 
development of rough surfaces, cellular patterns and, in 
particular under oblique incidence of the ion beam, the 
formation of more or less regular ripple patterns. For an 
overview on recent experimental results and theoretical 
approaches using continuum modeling we refer to Ma- 
keev et ali£ and Valbusa et alA 



During the last five years, however, spectacularly novel 
experimental results have been reported by Facsko et 
al^*2ii showing that GaSb and InSb semiconductor tar- 
gets eroded by Ar + ions under normal incidence can de- 
velop into a rather well ordered surface morphology with 
basically hexagonally arranged dot structures. Similar 
results have been subsequently reported by Gago et al& 
for Si targets under normal incidence and, more gener- 
ally, by Frost et al.V ! 11 for rotated InP, InSb and GaSb 
targets under oblique incidence (where as function of the 
inclination angle a variety of other patterns have also 
been observed). To theoretically explain the hexagonal 
ordering and taking advantage of earlier work by Elder et 
a j ( i2 4 i3 a i g0 an even earlier studyi^ on this subject), 
Facsko et have recently suggested a (stochastically 
extended) stabilized or damped Kuramoto-Sivashinsky 
(KS) equation for the dynamics of the surface height H . 



d t H = -v -aH + vV 2 H-D cS V 4 H+-(yH) 2 +ri (1) 

and might be considered as a useful continuum model 
for the ion-beam erosion under normal incidence since it 
can successfully reproduce the experimentally observed 
hexagonally ordered dot- type structures. The physical 
origin of the six terms on the rhs of Eq. is attributed 
to constant erosion velocity vq, dissipation, effective sur- 
face tension (Bradley-Harper mechanism 16 ), thermal and 
erosion induced diffusion, tilt dependent sputtering yield, 
and some stochasticity of the erosion process, respec- 
tively (cf. the original studjsi^ for details). Previous 
attempts^ii based on the (standard) KS equation, i.e. 
without the term —aH, have reproduced cellular pat- 
terns without convincing evidence of a regular hexagonal 
ordering. As noted by Facsko et al. at the end of their 
papeji^, there are two basic problems with the appli- 
cability of Eq.QJ to erosion processes since, as Eq.QJ 
stands, it violates the translation invariance in the ero- 
sion direction and, moreover, the physical meaning of the 
dissipation term —aH is not directly obvious. Facsko et 
;i1pU briefly argued that (i) translational invariance can 
be restored by assuming that —aH has to be replaced 
by —a(H — H) with H being the erosion depth averaged 
over the sample area and (ii) that the term —a(H — H) 
might be interpreted as an approximation of a newly sug- 
gested redeposition effect of the sputtered target parti- 
cles. Due to the nonlocal character of this term, how- 
ever, the systematic connection between this nonlocally 
extended damped KS equation and Eq.QJ is far from 
being obvious and the main reason for our paper. 

First, we argue on general grounds, by using balance 
considerations, symmetries and allowing for nonlocal de- 
pendencies, what the simplest functional form of the 
spatio-temporal dynamics of the morphology of ion sput- 
tered surfaces under normal incidence of the ions might 
be, if nonlocal terms as suggested by Facsko et al^ 
are taken into account. Second, we show that the re- 
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suiting non-local model equation can be mapped to the 
damped KS equation by the use of a temporally nonlocal 
transformation. By that we show how the stabilized KS 
equation suggested by Facsko et alii^ systematically fits 
into the theoretical framework of a continuum descrip- 
tion of eroded surfaces obeying all desired invariances. 
Third, we discuss in detail the subtle interplay between 
potentially nonlocal terms in the evolution equation for 
the eroding surfaces and the role of the KS-nonlinearity 
(V-ff) 2 . Finally, we also investigate some consequences 
on the general evolution dynamics in these equations. 

II. BALANCE EQUATION CONSIDERATIONS 

The starting point of our investigation is the general, 
stochastically extended balance equation for the spatio- 
temporal evolution of the eroded surface morphology 
H (x, t) measured perpendicularly to the initially flat tar- 
get surface with coordinates x = (x,y). Assuming that 
the particle density at the surface of the target material 
can be basically considered as being constanliiSii^, such 
a balance equation reads quite generally 20 

d t H = V- J H +F + ti. (2) 

and expresses the fact that temporal changes of the ero- 
sion depth H(pc,t) arise from two main contributions: 
(i) the detachment of target material leading to an in 
general inhomogenous and non-stationary, appropriately 
rescaled flux of eroded particles F (given by the number 
of eroded particles per time and surface area divided by 
the particle density) and (ii) the local rearrangements of 
the particles at the surface leading to a relaxational cur- 
rent Jh along the surface. Note that this balance only 
accounts for the target particles. The underlying driving 
of the erosion process, the flux of eroding ions /, enters 
indirectly into such a description: all terms on the right 
hand side of @ depend on / in a way that they vanish 
if the sputtering process is turned off, 1 = 0. Note that 
thermally activated processes that do not depend on / 
and might lead to a further smoothing of the surface af- 
ter the sputtering has been stopped are not taken into 
account. The spatio-temporal fluctuations 77 = rj (x, t) 
entering in J2J) mimic some stochasticity present in the 
erosion process and are usually assumed to be Gaussian 
white, i.e. having an average (??(x, t)) v = and a covari- 
ance (r)(x, t)rj(x' , f )}„ = 2DS(x - x.')6(t - t'). 

Considering periodic boundary conditions on an ap- 
propriately chosen, large enough sample area of the 
size L 2 and introducing the spatial average ~ (t) = 

(f/i 2 ) j L ^ 2 dx j^^ 2 dy... (being generally not equiva- 
lent to a stochastic average if taken on a finite area), 
the evolution of the mean erosion depth H(t) develops 
according to d t H(t) = F + rj which directly leads to 

H(t) - H(Q) = ( F dt' + ( rjdt' (3) 
Jo Jo 



with H(0) being the spatial average of the initial con- 
dition -ff(x, 0) (which is usually set to zero under the 
assumption of an initially flat target surface). Already 
at this stage it is clear that the mean erosion depth H(t) 
is, in general, not a linear function in t implying a con- 
stant erosion velocity, but a rather complicated function 
that integrates over the history of the stochasticity and 
of the eroded flux. 

In order to specify the admissible functional forms of 
the right hand side of Eq.(J2J), we reconsider the three 
fundamental symmetry requirements that are considered 
to be basic for the spatio-temporal evolution of surface 
morphologies? : (i) no dependence of J5J on the specific 
choice of the origin of time implying invariance of @ 
under translation in time, (ii) no dependence of @ on 
the specific choice of the origin of the coordinate system 
implying invariance of © under translation in the direc- 
tion perpendicular to the erosion direction, and (iii) no 
dependence of @ on the specific choice of the origin of 
the -ff-axis implying invariance of © under translation 
in growth direction. These symmetry requirements ex- 
clude any explicit dependence of V • Jh and F on the 
time i, the spatial position x and the erosion depth H, 
respectively. Following Facsko et al.'s argument, how- 
ever, an implicite functional dependence on H — H and 
on the spatial derivatives of H is still admissable. Con- 
sequently, the detachment contribution in J5| is quite 
generally given by F[VH,H — H] where [VH, ..} is the 
short hand notation for any derivative or combination of 
derivatives of H being compatible with the scalar char- 
acter of F or H. 

To proceed we (i) additionally apply invariance under 
rotation and reflection in the plane perpendicular to the 
erosion direction (which is suggested by the experimen- 
tally observed amorphorization of the target surface due 
to the erosion) , (ii) expand the height dependent term 
of F in a power series in H—H, i.e. fp = J2 n Pn{H—H) n , 
and (iii) perform a gradient expansion of the contribution 
f a of F that solely contains derivatives of H . The latter 
implies that the lowest order terms (up to forth order in 
V and second order in H) are given by V 2 H, (VH) 2 , 
V 4 H, V 2 (V 2 iT), (V 2 #) 2 ), and V ■[(VH)(V 2 H)}. Keep- 
ing only the lowest order contributions, the erosion flux 
is determined by 

F = F [l + MH -H) + ai V 2 H + a 2 (VH) 2 ] 

+ O (V 4 , H 2 , (H - ~H) 2 ) (4) 

with Fq being a constant and negative since this part 
of the flux is antiparallel to the erosion direction. The 
lowest order possible mixing term between Vi? and H — 
H that could appear in (0J is given by (H — H)V 2 H and 
will be omitted since it can be considered as a higher 
order contribution to the term proportional to oi\ or [3\ 
in JIJI. 

The functional form of the term representing the re- 
laxational currents at the surface, V • Jh in © remains 
to be specified. Following Facsko et ali 1 ^, we suppose 
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V • Jh = —D e ffV 4 H. This mimics the tendency of sur- 
face particles to reach energetically more favorable posi- 
tions with a positive curvature V 2 H > and, therefore, 
leads to a current Jh oc V(V 2 -ff). 

Combining all ingredients and renaming the entering 
coefficients b — /3iFo, a% — otiFo, a 2 — —D e g, and 
0-3 = CX2F0, yields a minimal functional form for the evo- 
lution of the erosion depth under the afore-mentioned 
restrictions that is given by 



d t H = ai V 2 H + a 2 V 4 H + a 3 (VH) 2 
+b(H -H) + F + r). 



(5) 



The functional form of Eq. © looks like the equation sug- 
gested by Facsko et ali^ that follows after their ad-hoc 
replacement of —aH by —a(H — H) in the damped KS 
equation. Note however, that the spatially constant part 
of the flux F, i.e. Fq, is, in general, not given by the mean 
constant erosion velocity since the latter is not constant 
during the course of the evolution. 

Next we transform Eq.© into a coordinate system 
that moves with the mean erosion depth H(t), 



h(x,t) = H(x,t) - H(t) 



(6) 



with h(x, t) being the local erosion profile that obviously 
fulfills h = and dfh = for all times t. The trans- 
formation © also guarantees that the invariance under 
translation in erosion direction H — > H + z with z be- 
ing an arbitrary constant (that only needs to hold in the 
coordinate system that is fixed in space) holds for any 
erosion profile h(x,t). As a consequence, the evolution 
dynamics for h(x, t) does not necessarily need to fulfill 
this requirement. 

Using Eq.© and the facts that dtH = dtH + dfh and 
that any terms consisting only of spatial derivatives of H 
have the same functional form in the comoving system 
when H is substituted by H + h, the evolution equation 
for the erosion process in the comoving frame reads ex- 
plicitly 

d t h = bh + a 1 \7 2 h + a 2 \7 i h + a 3 {\7h) 2 +F -d t H + r] (7) 

Eq. © still constitutes a stochastic integro-partial differ- 
ential equation by virtue of the nonlocal term dtH. Only 
if dtH = Fo or, equivalently, the mean erosion depth 
H (t) were moving with a constant speed for all times, 
Eq.Q would reduce to the damped stochastic KS equa- 
tion. In general, however, this cannot be invoked, in 
particular because of the presence of the KS-type non- 
linearity (Vh) 2 . By taking the spatial average of Eq.©, 
a 3 (V/i) 2 + Fo — d t H + rj = 0, one can directly connect 
the mean erosion depth H(t) to the dynamics of the local 
erosion profile /i(x, t), 



H(t) = F t 



dt' 



a 3 (Vh) 2 + r] 



(8) 



where H(0) = has been assumed. Consider first the 
deterministic part of i.e. with rj = 0. Since the 



integrand of the second term on the right hand side of 
Eq.© is positive for all times, the mean erosion depth 
H (t) systematically deviates from a linear time evolution 
given by F$t. For a 3 being negative (positive) the mean 
erosion depth is therefore retarded (advanced) in com- 
parison to Fot. Noteworthy, this fact is not a specific 
property triggered by the nonlocal term in Eq.JSJ), it is 
already present in the standard KS equation and from 
related studies therein, it is known that the second term 
of the right hand side of © significantly contributes to 
the time evolution of H(t). The additional impact of 
the term J Q dt'rj is roughly that of a superposed Wiener 
process since rj has been assumed to be a Gaussian and 
white. 



III. TRANSFORMATION TO A LOCAL 
EQUATION 

As explained in the previous section, it is not obvious 
how the nonlocal equations JSJ and J7J are related to the 
damped KS equation. Here we show in a more general 
context how the nonlocal term can be elimated by an ap- 
propriate transformation. Specifically our statement is as 
follows: There is a transformation h — > h that maps the 
general form of a nonlocal stochastic evolution equation 
given by 



d t h = bh + G[Vh] +F -d t H + r) 



(9) 



to a local evolution equation in the transformed variables 
h reading 

d t h = bh + G[Vh] + i] (10) 

where, quite generally, the functional G in the Eqs.© 
and l|l(J|) can contain any combination of derivatives of h 
or h, respectively, but no explicit dependence on h or h. 

As subsequently useful observation, we note that 
Ea. l|10(l is invariant under the transformation h — > h + 
z exp(6t) with z being an arbitrary constant. Obviously, 
@ contains the model equation J7J as special case if we 
set G[Vh] = ai V 2 h + a 2 V 4 h + a 3 (V/i) 2 . 

To the constructive derivation of the statement: In 
order to find the desired transformation we use the ansatz 

h -^h = h+^ + f(t) (11) 
b 

where the time-dependent function f(t) is so far arbitrary 
and will be subsequently determined. Since d t h = d t h — 
dtf and spatial derivatives of h transform without any 
change to spatial derivatives in h, insertion of (|ll|l into 
Eq-© yields 



d t h = bk+ G[Vh] + V + dtf -bf- d t H. 



(12) 



In order to arrive from Eq.JHJ at Ea. (jlOfl . one has to 
demand that the last three terms on the right hand side 
of Ea.ip.2p vanish at any time, 



dtf = bf + d t H. 



(13) 
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The linear nonautonomous ordinary differential equa- 
tion l|13|) can be solved for f(t) by standard means, e.g. 
by the method of variation of constants. Introducing 
f(t) = c(t) exp(bt) and solving the resulting equation 
d t c = exp(—bt)d t H yields 



/(*) = e' 



bt 



c(0) 



dt'e- bt 'd t >H(t') 



(14) 



In general, the integration constant c(0) in i|14|) can be 
arbitrarily chosen. This is a consequence of the afore- 
mentioned invariance of Ea. l|10|) with respect to h — > 
h + exp(bt)z (z arbitrary and constant) which implies 
that there is actually a continuous number of transfor- 
mations leading from Eq.@ to Eq. l|lUI) . A convenient 
choice, however, is to select as initial condition for h that 
/i(x, t) and h(x, t) initially coincide, i.e. /i(x, 0) = h(x, 0). 
Consequently, c(0) = —Fo/b holds and the transforma- 
tion reads 



h(x,t)=h(x,t)-^-(l-e bt ) 



dt'e- bt 'd t >H(t'^ 



(15) 

The transformation l|15|) that reduces the nonlocal equa- 
tion © in /i-system to a local equation (|10|) in /i-system 
has several specific properties, (i) It is a purely temporal 
(integral) transformation relating h(pc, t) in the physically 
meaningful coordinate system comoving to the mean evo- 
lution of the erosion depth with the evolution of h(x,t) 
in a temporally shifted system with no obvious physical 
significance. This shift s(t) = h(x,t) — h(x,t) being rep- 
resented by the last two terms in l|15|). takes over the 
nonlocal properties of Eq.@ and is therefore nonlocal 
by virtue of the last term in I] 15(1. Noteworthy, s(t) in- 
tegrates over the temporal history of the mean velocity 
of the erosion front via d t 'H(t'). (ii) The stochastic part 
in Eq.Q remains unchanged by the transformation, (iii) 
The transformation (|15fl is highly useful because it allows 
for the direct applicability of theoretical and numerical 
results obtained from Ea. (|10f> to Eq.@. Note, however, 
that for a correct interpretation in the physical meaning- 
ful H- or /i-system the full temporal information of the 
mean evolution of the erosion front H(t) needs to be sep- 
arately determined. If the dynamics of h(x, t) is known, 
this can be achieved by using 



H(t) = F t 



dt' 



G[Vh)]+rj(t' 



(16) 



(iv) The basic prerequisite for the transformation 1|15|) 
is that Eq.(|!5J) is linear in h. Extensions to nonlinear 
dependences on h in Eq.Q do not seem to be generally 
feasible. 

For the specific case under consideration, G[V/i] = 
aiV 2 h + ei2V 4 /i + a 3 (Wh) 2 , the transformation l|15fl can 
be further simplified. Using © it follows that 

h(x, t) = h(x, t)-a 3 f dt' e- i,(t '- t )(V/ip(t') (17) 



and, consequently, 

d t h = bh + ai V 2 h + a 2 V i h + a 3 (Vh) 2 + r) (18) 

which constitutes the damped KS equation in the h- 
system and possesses the obvious invariance under the 
transformation {h, a 3 } — > {— h, — a 3 }. Ea. lfTTjl shows the 
importance and the genuine interrelation of the KS-type 
nonlinearity for the non-triviality of the transformation. 
If 03 equals zero, then simply /i(x, i) = /i(x, i) follows. 

This argument can be generalized: Separating the 
functional G[V/i] in Eq.@ in terms that can be rewrit- 
ten as the divergence of a flux, Gp[V/i] = V • jp, 
and terms Gnf[V/i] that cannot, the corresponding 
transformation is determined by h(x,t) = h(x,t) — 
Jq dt' e~''^' _ ^GNF[V/i](t'). Consequently, any term be- 
ing nonlinear and not originating from a flux leads to 
analogous complications as the term (V/1) 2 . Only if 
Gnf[V/i] = 0, then again /i(x, t) = h(x, t) holds implying 
that ~H{t) = F t. 



IV. SOME FURTHER RESULTS 

A. Role of the nonlocal term 

The general idea behind model equations such as 
Eq.JSJ) is that any individual term entering into the sum 
on the right hand side has its individual physical sig- 
nificance. As argued by Facsko et al. 15 , the nonlocal 
term in JSJ) might be interpreted as a redepostion effect 
of eroded particles. To clarify the role of the nonlocal 
term in Eq.JSJ), b(H — H), we disregard, for the mo- 
ment, all terms in Eq.Q that depend on derivatives of 
H and the stochastic fluctuations 77. So, the equation 
dtH = Fq + b(H — H) remains to be solved. Obviously 
dtH = Fo holds and, consequently, (i) the mean depth 
evolution increases linearly with time for all times t ac- 
cording to H = Fot and (ii) H — H can be rewritten 
as H — Fot showing that this term, individually consid- 
ered, acts locally in the frame moving with F^t. The 
solution for the erosion profile H is then simply given 
by H (x, t) = H(x, 0) exp(fet) + Fot. For negative b, the 
impact of the term b(H — H) is just an exponential dimin- 
ishment of the initial target profile H(x, 0) in such a way 
that the overall shape and, in particular, the maxima, 
minima, and inclination points of -ff(x, t) — Fot remain 
at same spatial positions. Only the amplitude of -ff(x, t) 
decays in time. Consequently, it is not obvious whether 
the term b{H — H) can be interpreted as a redeposition 
effect that should also lead to a lateral variation of the 
erosion profile. 



B. Stability of a flat erosion front 

Taking advantage of the afore-mentioned arguments on 
the role of the nonlinearity in JSJ, ignoring the stochastic 
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fluctuations for the moment and linearizing Eq.JSJ, i.e. 
omitting the term (Vi?) 2 , yields H(t) = F t. Therefore, 
a flat erosion front given by -Hff(x, t) = Fot or, equiva- 
lently /iff(x, t) = solves the deterministic limit of the 
linearized version of Eq.JSJ) and determines its basic so- 
lution. Consequently, as far as the linear stability of the 
flat erosion front solution is concerned, the nonlocal term 
can be replaced by the local term F$t in the linearized 
Eq.0. Then, standard techniques can be invoked for the 
stability analysis of Hp-p: Using an ansatz for a perturba- 
tion of /iff( x j t) = of the form h oc exp[ik-x+(ji], yields 
a dispersion relation a = b — a\k 2 + a^fc 4 (fc 2 = k 2 ) for 
the growth rate <r(k) of perturbations with a wave vector 
k. Therefore, cr(fc) has its maximum at fc^ ax = ai/2a2 
with cr(fcmax) = b — a 2 /4a2 and the flat erosion front Hff 
is stable if b < a\/Aa,2- 



C. Time dependence of the mean erosion depth H 

A crucial point of our discussion is the fact that varia- 
ti ons of t he mean erosion velocity about Fo, dtH — Fq = 
a^(VH) 2 + r\, are generally non-zero and not even con- 
stant during the course of the erosion process. Here, we 
substantiate this by numerical simulations of Eq.Q for 
a representative set of parameter values that leads to 
hexagonal patterns and draw some further conclusions. 

In Fig. 1 we show the evolution of the ensemble or 
stochastically ave raged e volution of d t H — Fo given by 
d t {H) - Fo = a 3 ((VH) 2 ) since (rj} = {rj} = 0. For this 
quantity, the stochasticity necessarily present in individ- 
ual runs of the erosion process, as shown in Fig. 2, is 
leveled out. In both cases, the initial target profile has 
been kept fixed to a small Gaussian distribution of the 
initial amplitudes with a maximum amplitude of rj = 0.01 
about the perfectly flat state. 

The generic behavior as function of time consists of 
three parts: (i) For very small erosion times, dt(H) — Fo 
is very close to zero since here the nonlinearity in the evo- 
lution equation can be neglected and, therefore mainly 
the linear evoluti on of H contributes. Since dtH — Fo is 
proportional to asiVH) 2 this leads to a purely exponen- 
tial increase for short times as depicted in the insets of 
Fig. 1. (ii) With increasing erosion time, the amplitude 
H eventually reaches values where the nonlinear term in 
Eq.JSJ) becomes comparable in size to linear terms. In this 
crossover time range, dt(H) — Fo increases very rapidly 
up to a point where the subsequent increase drastically 
slows down, (iii) For longer times, where the nonlin- 
earity in Eq.(JSJl is fully developed, the KS-nonlinearity 
(ViJ) 2 mainly reduces the further increase of dt(H) — Fq 
almost to a constant. Our numerical simulations that 
went very far into the full nonlinear regime, however, 
do not indicate a full saturation, but a slight system- 
atic increase with time that depends almost linearly on 
t in the simulated time range. Similar behavior for in- 
dividual erosion processes (albeit slightly modified by its 
intrinsic stochasticity) can also be read off from Fig. 2. 



An important consequence of this dynamical behavior of 
Eq.JSJ is the fact that the underlying surface morphol- 
ogy -ff (x, t) does not saturate in a steady state unless 
dtH — Fo fully saturates into a constant value. To sub- 
stantiate this statement, we show in the insets of Fig. 
2 for two representative erosion times the corresponding 
morphology behavior that clearly indicates a non-steady 
evolution. 

Finally, we like to draw attention to a challenging point 
for further experiments. Due to the non-monotonic be- 
havior of dt (H) — Fo as function of the erosion time that 
is triggered by the KS-nonlinearity, a specific experimen- 
tal measurement of this quantity can lead to experimen- 
tal evidence for the existence and dominance of no-flux 
terms as discussed in section III. 



V. SUMMARY 

To conclude, we stress the most important points of 
our investigation: (i) The nonlocal, damped KS equa- 
tion suggested by Facsko et alp 1 ^ can be obtained as min- 
imal model from balance equation considerations if non- 
local terms compatible with the underlying invariances 
are allowed. Most importantly, such a nonlocal term be- 
ing proportional to H — H leads in combination with 
the KS nonlinearity to a non-monotonic time evolution 
of the spatially averaged erosion depth and, connected 
with that, to a surface morphology that does not gen- 
erally seem to approach a steady state (at least on the 
simulated time scales), (ii) The necessity for a non-local 
dependence arises from the fact that the KS nonlinearity 
cannot be rewritten in form of the divergence of a flux. 
If all terms entering into the evolution equation could 
be rewritten in form of such divergences, the non-local 
term could be replaced by a simple drift term H = Fot. 
(iii) By using a specific transformation presented in sec- 
tion III, the nonlocal, damped KS equation can be ex- 
actly recast in form of a standard (non-local) damped KS 
equation. This property greatly simplifies the analysis 
of the evolution equations. Finally, we have also briefly 
proposed a simple experimental measurement in order to 
substantiate the pronounced effect of the KS nonlinearity 
or at least related non-flux nonlinearities during the ero- 
sion process. Building on the results presented here, we 
will discuss, in a subsequent publication, the rich pat- 
tern forming structure of an appropriately generalized 
non-local anisotropic damped KS equation for the case 
of oblique incidence. 
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FIG. 1: Growth rate of the mean erosion depth calculated 
from Eq.0, averaged over 100 runs. Parameters are the same 
as in the paper by Facsko et al^&, except for a spatial step 
size of dx — 1 (mesh size 400 x 400,6 = —0.24, ai = — 1, 
d2 = — 1, <Xz = 0.0025, white noise with a maximum amplitude 
of X] = 0.01). Solid line: b = -0.22, dotted line: b = -0.24. 
The inset shows a semilogarithmic plot of a portion of the 
data where linear terms determine the growth. A linear fit in 
this regime yields ]n(d t (H)-F ) oc 0.018-t for b = -0.24 and 
ki(d t (H) -Fo) oc 0.057-i for b = -0.22. This compares well to 
the result which can be obtained from the linearized version 
of Eq. 0, i.e. 2a max = 0.02 for b = -0.24 and 2a max = 0.06 
for b = -0.22. 



FIG. 2: Growth rate of the mean erosion depth calculated 
from Eq.JSJ for five different noise sequences. Parameters are 
the same as in Fig. Qwith b — —0.24. Insets show a section 
of the surface morphology corresponding to the thick line in 
the figure at times t = 1 • 10 4 and t = 2 ■ 10 4 . 



